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Abstract 

In this paper we develop a procedure to compute the one-loop quantum correction to the kink 

masses in generic (l+l)-dimensional one-component scalar field theoretical models. The procedure 

uses the generalized zeta function regularization method helped by the Gilkey-de Witt asymptotic 

expansion of the heat function via Mellin's transform. We find a formula for the one-loop kink mass 

shift that depends only on the part of the energy density with no field derivatives, evaluated by means 

of a symbolic software algorithm that automates the computation. The improved algorithm with 

£\J respect to earlier work in this subject has been tested in the sine-Gordon and A(0)| models. The 

t-H quantum corrections of the sG-soliton and A(0 4 )2-kink masses have been estimated with a relative 

error of 0.00006% and 0.00007% respectively. Thereafter, the algorithm is applied to other models. In 

particular, an interesting one-parametric family of double sine-Gordon models interpolating between 

the ordinary sine-Gordon and a re-scaled sine-Gordon model is addressed. Another one-parametric 

O I family, in this case of </> 6 models, is analyzed. The main virtue of our procedure is its versatility: it 

can be applied to practically any type of relativistic scalar field models supporting kinks. 

PACS: 11.15.Kc; 11.27.+d; ll.10.Gh 
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1 Introduction 
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In this paper we shall address the broad topic of one-loop quantum corrections to kink masses in (1+1)- 
dimensional real scalar field theory models. Kinks are solitary waves that arise in non-linear models. This 
kind of solution resembles at the classical level an extended particle, i.e., they are localized, finite-energy 
objects. Some particles in nature, mostly baryons, seem to be extended, although they are sometimes 
approximated by point particles in theoretical calculations. In 1961 Skyrme pointed out that, in a certain 
extension of the non-linear sigma model now called the Skyrme model, there exist both 3D dispersive and 
solitary waves among their solutions [ ]. Because the model attempted to describe low-energy hadron 
phenomenology, and because solitary waves are formed from a heavy classical lump of energy, the idea 
is natural: upon quantization, dispersive waves become light mesons -pions- and solitary waves give rise 
to heavy baryons -protons, neutrons-. This bold idea prompted the task of investigating solitary waves 
in the quantum domain, a task mainly performed in the seventies, see e.g. [2, 3, 4]. 

The first theoretical studies in this issue were accomplished by Dashen, Hasslacher and Neveu by 
studying the quantum A(</>)2 and sine-Gordon kinks. In 1974 they succeeded in computing the one- 
loop correction to the classical mass of these solitary waves by developing the ^-expansion of these 
(l+l)-dimensional field theories. From the work of these authors in papers [2, 5] emerged the Dashen- 
Hasslacher-Neveu (DHN) formula 1 for the kink mass quantum correction. The DHN formula demands 
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To be precise, the formula for static kinks should be called the first DHN formula. There is a second DHN formula that 
calibrates the WKB correction to the Bohr-Sommerfeld levels of quantum breathers, see [34] . 



knowledge of the spectrum of the variational second-order (Hessian) differential operators valued at 
the kink/soliton solutions. In the s-G and A(0)| models these operators are Schrodinger operators of 
Poschl- Teller type, whose spectral problem is solvable. Moreover, the second-order fluctuation operators 
around the sG and A</> 4 kinks are respectively the first and the second of Posch- Teller Hamiltonians in the 
hierarchy of transparent - reflection coefficient equal to zero- models of this type. This provides sufficient 
spectral information to achieve exact computation of the one-loop quantum correction to the kink masses 
through the DHN formula: 1) A£ sG (0 x ) = -fia 2) AE^ifa) = ^ (^ - f)- 

Apart from these two models the DHN formula has little computational value. The reason is that the 
Hessian differential operator valued on the kink solution that emerges in other models has an unknown 
spectrum in general. Therefore, computation of the one-loop quantum correction to the kink mass in 
a generic model required the development of a new procedure; see the Reference [1 ] to find a summary 
of the work of our group on this issue during the last ten years. In fact, we started the new approach 
simultaneously with renewed interest about the problem of the quantization of classical lumps at the end 
of the last century. The new impetus came from the quantization of supersymmetric kinks, where the 
contributions associated with fluctuations of the bosonic and the fermionic fields balance each other, which 
avoids the need to identify the spectrum of the Hessian operators involved in this problem. However, 
new subtleties concerning the effect of the boundary conditions come into play. Several groups at Stony 
Brook/Wien, [8, 9], Minnesota [10] and MIT [11, 12] addressed mixed issues in the problem by studying 
the impact made by using different types of boundary conditions - PBC, Dirichlet, Robin-, regularization 
methods -energy cutoff, mode number cutoff, high-derivatives-, and/or performed phase shift analysis, 
in connection with possible modifications due to the quantum effects of the central charge of the SUSY 
algebra [13, 16]. The Stony Brook/Wien group formed by Rebhan, van Nieuwenhuizen, and Wimmer, 
together with Goldhaber, investigated the computations of mass shifts induced by one-loop fluctuations 
on supersymmetric kinks [17, 18, 19]. The extreme elusiveness of this issue did not prevent those authors 
from identifying the old DHN formula as being based on a regularization method that sets a cutoff in 
the number of fluctuation modes to be counted, rather than the conventional energy cutoff. Another 
group from Minnesota University addressed the same problem by using high-derivative regularization, 
with SUSY being preserved by boundary conditions to find similar results. Phase-shift analysis by an 
MIT group -Jaffe, Graham, and collaborators- also led to some advances, in this case in a purely bosonic 
setting [11]. 

Our goal in this paper is threefold. First, we construct a better structural derivation of the DHN 
formula and provide the one-loop kink mass shifts of the sG- and A</> 4 -kinks after a regularization by 
means of a cutoff in the number of modes. The same results are achieved by a treatment of the DHN 
formula based on the heat kernel/zeta function regularization of ultraviolet divergences. 

Second, the method of zeta function regularization was invented by Dowker and Critchley and, inde- 
pendently by Hawking, circa 1976 [ ]. The partition function of Euclidean quantum field theories is a 
functional integral that, up to one-loop order in the ^-expansion, is the inverse of the square root of the 
determinant of a differential operator of Laplace type times the exponential of the classical Euclidean 
action over h. In our framework, the pertinent differential operators are those ruling small quantum 
fluctuations in the classical kink backgrounds. Generically, the spectral information in these situations is 
grossly insufficient for identifying the generalized zeta function in terms of known spectral functions. For- 
tunately, in the mid sixties de Witt [2 ! ] had already proposed the use of the high-temperature expansion 
of the kernel of the generalized heat equation related to the differential operator of Laplace type to unveil 
the meromorphic structure of the generalized zeta function. To achieve this goal, one takes advantage 
of the link between generalized heat and zeta functions via Mellin transforms, such that the residua at 
the poles of the generalized zeta function are proportional to the Seeley coefficients of the heat kernel 
expansion [22, 23, 24, 25, 26, 28, 29]. We shall use all this theoretical machinery in order to construct a 
formula for the one-loop quantum correction to the kink mass in a generic (l+l)-dimensional real scalar 



field theory as a truncated series in the Seeley coefficients of the Hessian heat function. 

Third, the effectiveness of this method relies on the computation of the Seeley coefficients. This is 
not an easy task and requires the solution of difficult recurrence relations between the Seeley densities 
and integration of these densities over the whole real line. In this paper we shall improve our symbolic 
method implemented on a PC computer in a Mathematica environment in two ways: 1) We shall adapt 
in an optimum way the range of the inverse temperature in the integral Mellin transform to the number 
of Seeley coefficients kept in the truncated formula. 2) We shall use the first-order equations to skip the 
integration over the whole spatial line traded by integration in the field space between the two vacua 
connected by the kink. These two advances offers not only great computational advantage but also 
higher precision. The error in the computation of the one-loop sG-kink is diminished from 6.0% to the 
0.00006% and for the A0 4 -kink the reduction is from 0.07% to the 0.00007% !!! But this is not all: several 
other models for which the spectral information is not available can now be treated and the one-loop 
mass shift calculated. Among them the double sine-Gordon model [37, 38], the Razavy potential [39], 
the BGLM model [40] and others. Especially interesting is the application of the method to a family of 
trigonometric models interpolating between the sine-Gordon and the re-scaled sine-Gordon models. A 
family of deformations of the <f> e model [41, 42] will be also analyzed. 

The organization of this paper is as follows. In Section §. 2 and §. 3 we shall review the standard 
lore about kink topological defects and their one-loop fluctuations. This review will allow us to fix 
the notation and to remark on some important features of this kind of solutions. Section §. 4 will be 
devoted to the DHN formula. A cutoff in the number of modes and, alternatively, the generalized or 
spectral zeta function procedure will be used in the regularization of intermediate divergent quantities. 
We shall apply this formula to the (ft 4 and sine-Gordon model examples. In section §. 5 the Seeley-Gilkey- 
deWitt-Abramidi heat kernel expansion will be introduced to express the generalized zeta function as an 
asymptotic series. We stress that the computation of the Seeley coefficients only requires the form of the 
potential. Again, we address the Xcfr and sine-Gordon model examples in order to check the accuracy 
of these formulas. Section §. 6 is devoted to obtaining the asymptotic formula to be summarized in a 
symbolic algorithm, whereas in Section §. 7 we construct an efficient program for that algorithm. We 
test the algorithm on the c/> 4 and sine-Gordon models, gaining a feedback about the method. Finally, in 
Section §.8 we shall apply the program to several interesting scalar field theory models and families of 
models. 

2 Classical field theoretical models 

The action governing the dynamics in our (l+l)-dimensional relativistic one-scalar field theoretical mod- 
els is of the form: 

Here, ipiy^) '■ R ' — > R is a real scalar field; i.e., a continuous map from the 1 + 1-dimensional Minkowski 
space-time to the field of the real numbers. y° = r and y 1 = y are local coordinates in R ' , which 
is equipped with a metric tensor g^ = diag(l, -l),/j,,v = 0,1 such that y^ = g^y^yu, gf^-g§r = 

We shall work in a system of units where the speed of light is set to one, c = 1, but we shall keep 
the Planck constant h explicit because we shall search for one- loop corrections, proportional to H, to the 
classical kink masses. In this system, the physical dimensions are: 

[h] = [S] = ML , [y„} = L , M=M3L5 , [U] = ML' 1 . 

The models that we shall consider are distinguished by different choices of the part of the potential 
energy density which is independent of the field spatial derivatives: U[tp(y^)]. In all of them there will 



be two special parameters, mj and 7^, to be determined in each case, carrying the physical dimensions 
dimensions: [mj = L^ 1 and [7^] = M~2L~5. We define the non-dimensional coordinates, fields and 
potential in terms of these parameters: 

72 - 
x^ = m d y^ , x = t , xi = x , 6 = jaip , U(6) = —oU(ip) 

m d 

The action is also proportional to a non-dimensional action 

m = £s W = i // d*v Q * . £ - t/wx-)]) in 

and the field equation in this non-dimensional setup is the PDE: 

d 2 6 d 2 6 _ dU 
dt 2 dx 2 86 



(2) 



Unless U(6) is a quadratic polynomial in the field, the PDE (2) is non-linear. Moreover, we shall assume 
that U{6) is a non- negative twice-differentiable function of 6: U(6) E C 2 (IR) and U(6) > for E R. 

2.1 The configuration space and the vacuum orbit 

The configuration space C of the system is the set of field configurations -the continuous maps from R to 
R for fixed time t = to- for which the potential energy functional 
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is finite. Thus, C = {6{t ,x) E MapefR 1 ,R) / E[<f>] < +00}. Let 

M = {6^ J U{6®) = 0} (4) 

be the set of zeroes of U, which we assume to be a finite or infinite discrete set. The initial conditions 
6(0, x) = 6^\ 6(0, x) = provide the static and homogeneous solutions to the PDE (2). When i > 1, 
a process of spontaneous symmetry breaking is triggered by the choice of vacuum in the quantization 
of these classical systems. In our models there will be a discrete symmetry group G and the set Ai 
is the union of one or several orbits of the action of G on any of the vacua: Ai = \_\ A -k-, where 
Ga,A = 0, 1, • • • ,r, is the little group in the A orbit. The vacuum moduli space is in turn the set of 
orbits: M=^ = {(/> (0) , (1) , • • • , 6^}. 

Each field configuration in C must comply with the asymptotic conditions 

lim 6(t, x)eM , lim M^2 = . (5) 

x— >-±oo x— >±oo OX 

in order to guarantee the finiteness of the energy (3). Thus, the configuration space is the union of 
topologically disconnected sectors 

— i,j ij 
where Cu stands for the set of the configurations that asymptotically connects the element 6^ l > with the 
element 6^> in M, i.e., we assume that: 

lim 6(t, x) = 6 {i) E M , lim 6(t, x) = 6 {j) E M 

x— >— 00 x— »+oo 

Because temporal evolution is continuous (a homotopy transformation) , the asymptotic conditions do not 
change in time and the sectors Cij are topologically disconnected. Note that the elements of Ai belong 
to the Ca sectors. 



2.2 Kinks and solitary waves 

In the Ca+i sectors there cannot exist static homogeneous solutions. In these sectors there are, however, 
spatially dependent static solutions of (2) of solitary wave type: non dispersive non-linear waves, see e.g. 
[30, 31, 32]. For static field configurations the PDE (2) reduces to the ODE: 

dx 2 d<\> w 

The ODE (6) can be seen as the equation of motion of a mechanical system: understand <f> as the particle 
coordinate, x as the particle time, and — U as the particle potential energy. The second order ODE (6) 
admits the first-integral 

which is the particle energy in the analogous mechanical system. The asymptotic conditions (5) guar- 
anteeing finite field theoretical energy -particle action- force zero energy 1 = and the solitary wave 
solutions satisfy the first-order ODE 

g = ±^?M . (T) 

Thus, the solitary waves are in one-to-one correspondence with the separatrix trajectories between the 
bound and unbound motions of the mechanical system connecting contiguous unstable equilibrium points. 
In fact, the zero mechanical energy Hamilton characteristic function 

W(4>) =± ! d<f> ^2U{4>) 

provides the particle action -henceforth the solitary wave or kink energy - as a topological quantity: 

E{(j) K )= W{<f>®) - W{4> {i+1) ) . 

Moreover, a Lorentz transformation sends the static solution 4>k{x) -in the kink center of mass- to another 
time-dependent solution cf>K(t,x) = 4>k{ /i-^ f)) showing the kink as a traveling wave. 

2.3 The prototypes: sine-Gordon and </> 4 kinks 

We now collect very well known facts about these structures in the prototypical sine-Gordon and </> 
models for the sake of completeness. 

2.3.1 The sine-Gordon model 

In this famous model, the dimensionless potential in (1) is the trigonometric function 

U(4>) = 1 - cos 

whereas the special parameters are rr? d = m 2 and 7J = — j. Therefore, M = {<ft( n > = 2irn,n G Z}. The 
symmetry group is G = 7h<i x Z = B^, the infinite dihedral group encompassing the Z2 group generated 
by reflection in the fields <p ~^ ~ ( t ) times field translations in 2irn: cj) — )■ <j> + 2-7rn , n E Z. Thus, the set of 
zeroes of U is the orbit of the symmetry group, Ai = <j-, where in this case Go =1>i- In our context, 
the vacuum moduli space is the quotient space M = ^$ and it is a single point because the action of G 



on M. is transitive, meaning that all the elements in Ai are equivalent J\f = {<ft(°'}. Finally, the soliton 
or kink solution in the topological sector C nn ±i of (7) is: 

4> K {x) = ±4 arctan e x ~ x ° + 2im . 

The zero mechanical energy Hamilton characteristic function W{4>) = ±4cos| leads us to the kink 
classical energy: E(4>k) = 4 cos^n cos- 

2.3.2 The A0 4 model. 

In this model the dimensionless potential in (1) is a polynomial: 

U^) = ^ 2 -l) 2 

and the special parameters are m 2 , = *5- and 7J = -\. The set of zeroes of the potential is Ai = {(jr- 1 ' = 
(— l) l ,i = 0,1}. The symmetry group is G ~ {(f) — >• (f),(f> — >■ —(/>} = Z2, the discrete group generated 
by reflection in the fields. The set of zeroes of U is the orbit of the symmetry group <=r- = Ai where 
Go = I is the identity. The "moduli space of vacua" is again a one-element set, and the quotient space 
is: J\f = 4r = {4>(°>}. The kink solutions of (7) are now: 

4>k(%) = ±tanh(x — xq) 

The zero mechanical energy Hamilton characteristic function is W{4>) = ± ( ^ — <f>) . Therefore, the kink 
classical energy is: E(4>k) = II — 2 1 - 
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3 One-loop quantum fluctuations 

3.1 Vacuum fluctuations 

The small (quadratic) fluctuations around any of the equivalent constant solutions (j)(t, x) = 4>^' + 5(j)(t, x) 
satisfy the linearized field equations: 



JL_^ + ^ (( , I) + 0[W)2] = „ , || 



= v> . (8) 

0(0 



The solutions of (8) of the form 5<fik(t,x) = e tu ^ ' l fk{x) are the vacuum normal modes of fluctuation 
coming from the eigenfunctions of the Kq = — 4-$ + v 2 differential operator: 

K f k (x) = v 2 (k)f k (x) , u 2 (k) = k 2 + v 2 , f k {x) = e ikx , k£R . 



In order to tame the difficulties of dealing with a continuous spectrum, we place the system in a normal- 
ization interval of very large, but finite, non-dimensional length, I = m^L, and choose periodic boundary 
conditions on the plane waves: /&(— |) = /fc(|)- The spectrum of Kq becomes discrete and the eigen- 
functions are monochromatic plane waves: 

/t„W = -j e ik " x , k n = ^-n , neZ , v 2 (k n ) = -^-n 2 + v 2 . 

Therefore, the general solution of (8) with PBC is a Fourier series, the more general linear combination of 
these normal modes. In quantum theory, Fourier coefficients become creation and annihilation operators 
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of the fundamental mesons moving in the vacuum, and the energy of the vacuum state, when all these 
fluctuation states are unoccupied, is the one-loop shift to the vacuum energy: 

AE [^] = ^ A EotfO] = ^- ( lim f; u(k n )) . 



'YJ 2 \ N— >oo 



-N 



Putting the system in an interval with periodic boundary conditions, the vacuum spectrum becomes 
discrete and setting an integer number N to be large, but finite, a cutoff in the energy of the fluctuation 
modes is selected. At the end of the day, however, we must come back to the I — > oo and N — >• oo limits. 
Bearing in mind the spectral density of the fluctuation modes po(k) = 4? = ^-, we have the following 
expression for the vacuum energy at this double limit: 

AW W ] = ^ (/" dkpo(kHk) + |) = M (±1" dk V&^+ |) . 

We now remark an obscure and subtle point. At I = oo, a half-bound state with energy v^ = v and 
eigenfunction constant emerges. This state is different from the Vk=o = v state, the threshold of the 
continuous spectrum. The name is suggested by the fact that a factor multiplies the eigenvalue of the 
half-bound state in its contribution to the vacuum energy by one-half because of the one-dimensional 
Levinson theorem. To see a proof of this delicate point, see the APPENDIX in [33] which is worked 
following the seminal paper of Barton [14]. We also notice that in the I — > oo limit one goes beyond 
the Hilbert space L 2 (M) to its completion L 2 (M) by admitting monochromatic plane waves that are not 
square-summable. 

3.2 Kink fluctuations 

The small fluctuations over any of the equivalent kink solutions <j>{t,x) = <Pk{x) + 8(j)(t,x) satisfy the 
linearized field equations: 



r) 2 8 2 \ f) 2 TI 

w -^ + v* + V(x))5cf>(t,x) + O[(Sci>f]=0 , w 



= v 2 + V(x) . (9) 

(j) K (x) 



The solutions of (9), of the form S(pk(t,x) = e luJ ^ t f q (x), are the normal modes of kink fluctuations 
coming from the eigenfunctions of the K = —4^ + v 2 + V(x) differential operator: 

lim V(x) = , Kf q (x)=u J \q)f q {x) , oo 2 (q) = q 2 + v 2 . 

x— >±oo 

Because of conditions (5), the kink solutions asymptotically tend to the vacuum solutions when x —> ±oo. 
For this reason V(x) tends to zero both at x — > +oo and at x — > — oo. This is true in every model for 

a property that always happens if the vacuum moduli space is a set of 
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= y 2 = a2u 
4,(i) ' 9 <1> 2 

one element. In this case 2 , 



(i+l) ' 



lim K = K 

x— »±oo 



Unlike Kq, however, which is a Helmoltz operator, K is an Schrodinger operator. The kink normal 
modes are not plane waves but some dispersive wave functions distorted by the kink. Bound states and 
half-bound states can also exist. Thus, there are three types of eigenfunctions of K at the I = oo limit 
that are characterized as follows: 



2 In fact, lim.-,.,.,. K = K i+1 = -£ J+v ^ +1 , lim^oo K = K t = -fa +vf if %£\ = v\ / |^ = v? + i- This 



1 0(0 °1> |0(* + 1) 



happens when the kink interpolates between different points of the vacuum moduli space: cf> { >'' and (f) l+1 . 



1. Scattering states. The asymptotic behaviour is: 

f q (x) ~ exp[iqx + -S(q)] , q£R , 

where S(q) is the total phase shift induced by V(x). Periodic boundary conditions require that 

q n l + 5{q n ) = 2irn 

Therefore, q n = k n — \5{q n ) and the spectral density over the kink background is: p(k) = 5^ + 
1 dS(k) 
2-k dk ■ 

2. Bound states: / W1 (x), /^ 2 (x), • • • , fw b (x) with discrete positive eigenvalues cof = < wf < w| < 
• • • < oj 2 < v 2 below the threshold of the continuous spectrum v 2 and classified by the number of 
nodes b. The asymptotic behaviour is: 

lim /^(x)oce^-^ , lim /„.(*) a e^ v ~ Ui > , i = l,2,---,b . 

X— ►— 00 X— > + oo 

These wave functions are strictly L 2 (M). 

3. A half-bound state. If the value of the last eigenvalue is u 2 = v 2 then this eigenvalue involves 
the existence of a half-bound state, which is buried in the continuous spectrum threshold. This 
happens when the reflection scattering amplitude is zero. 

In particular, the lowest bound state in the kink sector is always a bound state. From ODE (6) one 
immediately sees that: 



d 2 d6 K d 2 U 



dx 2 dx 



>K 



. dx 

<PK 



There is always a zero energy eigenfunction of the K operator in the kink sector: f\(x) = -%£■ , u>i = 0. 
This eigenfunction is no more than the Goldstone boson due to the spontaneous symmetry breaking by 
the kink of the invariance of the system with respect to spatial translations: x —> x+a. Moreover, because 
4>k(x) is a monotonic increasing (or decreasing) function of x, the translational zero mode f\{x) = -^ 
has no nodes. According to Sturm-Liouville theory /i(x) is thus the lowest energy eigenfunction - the 
ground state - of the non-negative Schrodinger operator K. 

We write the general solution of (9) in terms of these kink normal modes, assuming a normalization 
interval of length I and PBC 3 : 



5(f>(t,x) 



00 1 1 

n^oo VMVn) V ' J 

We have not included in this expansion the zero mode because this mode does not enter at one-loop order 
(its existence is due to the spontaneous translational symmetry breaking). Also, the possible half-bound 
state arises only at the I = 00 limit when the energy of the highest bound state is exactly the energy of 
the scattering threshold. In that case, the Levinson theorem also sets a ^-weight to this state; see [33]. In 
quantum theory, the coefficients of the expansion become creation and annihilation operators of mesons 



3 Quantum field theories that are based on one-particle states containing some bound state are non-unitary. This is, 
however, the situation that we encounter in the kink sectors, see e.g. [35] 



moving either on the kink classical background -the scattering states-, or trapped by the kink: the bound 
states. The shift in the kink classical energy due to these unoccupied one- loop kink fluctuations is: 



AE ^ K ] = -± A E [<t> K ] = -^ [ Jim^ J^ w{q n ) + £ ^ 

-N j=2 



N b 



Sending the length I of the normalization interval and the cutoff in the energy iV to infinity, this formula 
becomes: 



fry 2 / Z" 00 ^ 

AE [<f) K ] = -£■ / dkp(k)u(k) + J2 "j + s b^b 



H ' /x 



(/I d Ki + ^ ( "O v/ ^^ + g a;i + S6a;6 ) ' 



2 
where s;, = 1 if uj\, < v and s;, = \ if ^ = v. 

4 Kink mass quantum correction. The DHN formula 

It is clear that both AE{(/)( 1 '] and AE[<Pk] are divergent quantities (series or integrals) because there is an 
infinite number of fluctuation modes. In the mid-seventies, Dashen, Hasslacher, and Neveu implemented 
a regularization and renormalization procedure in two steps. 

4.1 Zero-point energy renormalization: mode by mode subtraction 

In the first step, a zero point energy renormalization is performed by means of a mode-by-mode substrac- 
tion. This kink mode-by-mode energy subtraction can be written symbolically as: 



A£iN = ^Tr L2 \K*-K*j = -f^X, b>i-"i) • ( 10 ) 

where the index j = 1,2, • • • , N labels the eigenvalues ordered in increasing energy of both operators, 
when the spectra are discrete because some boundary conditions are chosen at the frontier of some finite 
interval of length I. In the Appendix of [ ], however, we showed, using the one dimensional Levinson 
theorem contextualized in the kink and vacuum fluctuation spectra described above, that the subtle 
I = oo limit produces the formula: 

A * ) = T fv ■^■V / ^ + ^+^(V(x)) + |J Wi + Sb a, 6 -Hj , (11) 

where (V(x)) = j_ dxV(x). The two last terms in (10) collect the contribution of the possible half- 
bound state of K and the ever present half-bound state of Kq. In the sG and (j) 4 models, the net effect 
is null but we shall consider other kinks such that there are no half-bound states in the spectrum of K, 
thus generalizing the DHN formula to the form shown in (12). We remark that, in our derivation of 
this formula, again following [' !], we relied on even and odd states complying with Dirichlet boundary 
conditions. Although not usual in QFT, this approach is convenient for two reasons: 1) There is no 
degeneracy between eigenstates with opposite wave number, which helps to distinguish the half-bound 



state from the continuous threshold. 2) The even state of energy v 2 has wave vector |, which gives rise 
(rather indirectly) to the ^-weight. 

In the standard literature on the subject, the divergent integral is regularized by means of a cutoff: 

ae 1 ( ( /> k )[a] = M | j^ d ± . ^ . y^r^ + ±(v(x)) + £ «; + ^ - w 

such that AEi(<j)jc) = limA->oo AEi {(J>k) [A] . We say that the kink Casimir energy is regularized by 
means of a cutoff A in the number of modes. 

4.2 Mass renormalization 

AEi(4>k) in (11) is still divergent because in models with interactions there are more than the vacuum 
energy ultraviolet divergences. In (1 + l)-dimensional scalar field theory normal ordering tames all the 
uv divergences. The contraction of two field operators at the same point is equal to the normal ordered 
product plus a divergent term coming from the one-loop graph in Figure 1 below: 



«^V):w , ^-kx^^rjt^ 



v' 2 



where we have used periodic boundary conditions /(— |) = /(4) -henceforth plane waves as one-particle 
states- to place the system in a normalization interval of non dimensional length I as usual in QFT 4 . 




Figure 1: Self-energy graph due to the valence-four vertex with factor < tM-((j)^'). 
The Hamiltonian density operator 

#[&*")] = \ U\ x n + 1^(^) ) + uifar)) , 



normal ordered in the standard form (Wick theorem) 
;%■=%+; ( 1-exp 



U(4,{x")) ■= W-^A. 2 : |^W(x")) : +0 H^iSv' : |£(0 : 



^> 2 4 



shows the necessary one- loop counter term 5 . Therefore, the self-energy counter-term contributes to the 
one-loop kink mass shift: 



f^T d X ((0;K\: d ^):\K;0)-(0;V\:^ 



oo 



AE 2 {4> K ) = -^6v 2 / dx I(O ] K\:—(0):\K;O)-(O;V\:— T (4>):\V;O} 



4 Henceforth, we shall use PBC. At the I — oo limit, however, we shall incorporate the contributions of the half-bound 
states derived in the ultra-long distance regime of the spectrum with Dirichlet boundary conditions. 

5 In 2D scalar field theory there are only ultraviolet divergences in closed-loop diagrams with only one internal line. 
Normal ordering is enough to control ultraviolet divergences. 
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Here, |0;.ff) and |0; V) are respectively the ground states - no mesons at all - in the kink and constant 
vacuum backgrounds. These states are coherent states and the action of normal ordered operators in this 
type of quantum states is easily calculated: 

(0; K\ : W (j>): \K; 0) - (0; V\ : ^) : \V; 0) = ^(<M*)) - W ^) = V(x) . 
Consequently, 

and the DHN formula for the one-loop kink mass shift reads: 

AE^k) AE 1 {cj> K ) + AE 2 {4> k ) l£i 1 v 



El v 

Uj + -s b cv b - - + (12) 



fr>2 H 2 U" J ' 2 

1 f°° „ 8S(k) r- 1 ^ / 1 /■« dfe 



which has been obtained from the regularized version when the cutoff goes to infinity: 

AE(<p K ) AE(4> K )[A] 1^ 1 v 

-HJT = A" w d = 2 g^ + 2^- 4 + (13) 



+ — lim / dk^^Vk 2 + v 2 + ^- (V(x)) (l-- lim / 

2irA^ooJ dk 4vr x v ;/ V 2 A->oc ,/_ 



2 A->oo 7_ A ^/A; 2 + u 2 

4.2.1 The prototypes: the sine-Gordon and A0 4 kinks. 

In the sine-Gordon model we have: 

d 2 d 2 

K = -— + 1 , K = -—, 7 + l-2sech 2 x , (F(x)) = -4 
aar aar 

SpecKo = {vi = l}iU{is 2 {k) = k 2 + l} keR 

SpecK = {oJi = 0}U{uj 2 = l}i\j{u> 2 (q)=q 2 + l} qm 

I ., . n 1 J 1 d%) / 1 1 
/°o(«0 = JT" > % =2arctan- , p(q = — + ^— = 7^ 1 , 2 

27T g 2tt 2tt dq 2i 7rl + q z 

Discarding the zero mode in the spectrum of K the eigenvalues of K$ and -fC are the same: one half-bound 
state and a continuous spectrum starting at the threshold lo 2 = 1. Because of the phase shifts induced 
by the kink potential well the spectral densities are different. A mode- number cutoff regulator applied 
to the DHN formula provides the exact result: 



Note that a cutoff in the energy would leave — > 2 = as the one-loop kink mass shift. 
For X(p 4 kinks things are sligthly more involved: 

d 2 d 2 

K = --j-^+4 , K = -- T +4-6sech 2 x , (V(x)} = -12 
dx z dx z 



11 



Spec Kq = 


- {HL = 


= 4}iU{^)=fc 2 + 4} M 




Spec K = 


= {wi = 


= 0} U {a| = 3} U {w§ = 4}i U {w 2 (g) = q 2 + 4} geR 




Po{k) = 


z 

2^ 


.. . 3q I \ d6(q) 
, 5(q)= 2arctan 2 _ (?2 , p(g) = ^ + ^ ^ = 


Z 

" 27 



vr(l+g 2 )(4 + (7 2 ) 

In this system, besides the zero mode there is another bound state in the spectrum of K with eigenvalue: 
u>2 = 3. All the other eigenvalues of Kq and K are the same: one half-bound state, to 2 = 4, and a 
continuous spectrum starting at the threshold: oo 2 (q) = q 2 + v 2 . The phase shifts induced by the kink 
potential well cause different spectral densities as in the sine-Gordon model. The mode-number cutoff 
regulator applied in the DHN formula also provides the exact result: 

AE( M v/3 3 f* 2 + k* 3/ M 1 



^TJ 2 7T A^oo 7 v / (l + A: 2 )(4 + /c 2 ) vr V A^°° y y/tf + 4 

\/3 3 .. / . , ,A. 1 V3A A 3 / A 

= lim arcsmhl— ) -\ — arr.tani 



..... . arcsinhf — ) -\ — arr.tanf - — ) ( 1 — lim arcsinhf- 

2 vrA-^oo^ K 2' ^3 V4 + A 2 J vr \ A^oo v 2 

\/3 3 1 13 

2 ir V3~ 2^/3 vr 

A cutoff in the energy provides the answer: — r % K = ^75- 

5 Zeta function regularization procedure 

The problem with the DHN formula is that one needs complete spectral information about the K oper- 
ators governing the kink fluctuations and, at present, it is available only for the sine-Gordon and \{(f)2 
kinks. Instead, we shall rely on the spectral AT-heat kernel expansion for computing one-loop kink mass 
shifts in other interesting models where there is not enough spectral information about the kink fluctua- 
tions. However, before this we shall regularize the DHN formula by the zeta function procedure to test 
the conceptual nature of this method on sine-Gordon and A0 kinks. 

The traces and determinants of powers of elliptic operators are only defined by means of a process of 
analytic continuation that mimics the definition of the Riemann zeta function as a meromorphic function, 
provideng formal meaning to strictly divergent series in some regions (Res < i) of the complex plane. 
By replacing natural numbers by eigenvalues (hopefully forming a discrete spectrum), generalized zeta 
functions associated with differential operators can be defined. If we assume that A G „4[L 2 (M)] is an 
operator with a positive discrete spectrum in the form of Spec A = {a; 2 G E : Aip n = uj 2 ip n }, the 
generalized zeta function associated with the operator A is a meromorphic function in the complex plane 
defined as 

( A (s) = TrA-° = Y,Ass ' seC 

There is a general theory of elliptic pseudo-differential operators that characterizes the conditions under 
which the generalized zeta functions are meromorphic functions, and the values apart from the poles 
of the zeta function and the derivatives of zeta, are taken as "regularized" definitions of traces and 
logarithms of determinants of (complex powers of pseudo-) differential operators [22, 23, 24, 25, 26]. 

Direct computation of the generalized zeta function is difficult in general. In order to identify the 
poles of (a(s), we introduce a non-dimensional fictitious temperature (3 = jp?jk 6 and define the A-heat 



Note that, in our system of units, the Boltzmann constant ks and fictitious temperature T respectively have dimensions 
of ML andL" 1 . 
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function: 

h A {P) = Tr e~P A = Y^ e _/3 ^ (14) 

n=0 

which has better convergence properties. Note that for operators with a positive spectrum the low- 
temperature asymptotic condition 

lim h A (l3) = 

holds. The Ca( s ) an d the h A ((3) functions are related through the Mellin transform: 

1 



Us) 



T(s) 



dpp s -'h A {p) 



(15) 



5.1 Vacuum spectral zeta and heat functions 

We shall now consider the L 2 trace of the Kq operator using periodic boundary conditions on the normal- 
ization interval. In the previous Section, we mentioned that Dirichlet boundary conditions are convenient 
to disentangle the half-bound state that arises at the I = oo limit from the threshold of the continuous 
spectrum through the one-dimensional Levinson theorem. We now know how to correct the kink Casimir 
energy to take this subtle effect into account. 

The energy of the free Hamiltonian evaluated at the vacuum state is: 



A%«) = (0;V\mV;0) = ^Tr L2(§1) 4 = ^J2 



niz- 



4f 



-n + v 



hm d 1 N 

2 W 2' 



Here, H 2 is the free Hamiltonian; i.e., the quadratic Hamiltonian in the fluctuations around the vacuum, 
\V;0) is the vacuum state where all the fluctuation states are unoccupied, and (k (— \) is the general- 
ized/spectral zeta function of the differential operator Kq at the s = — t> point in the s-complex plane. 
In fact, one considers in general the trace of Kq S 



(k (s) = Tr L2(§i) K s = ^2 



1 



E[s,v 2 



P 



s G 



r47r 2 „2 i „,2 

which is a convergent series only for Re(s) > ^ but can be analytically continued to the whole s-complex 
plane to find the meromorphic Epstein function, E[s, u 2 |-S-]. We also define the heat function: 



h K {P) = Tr L 2 ( gi)( 



-f3K 



-v 2 H 



E 



-P^fn 2 



-We 



nei 



(0|i£/J) 



where 



a 
b 



ZT 



E„ e z e 2lTi[in+a)iz+b)+{J ^ T)] (a, 6 = 0, \) in the form of a series that is a Riemman 



Theta constant: a Riemman Theta function at z = 0. The spectral K^-zetSi and i^o-heat functions are 
related via the Mellin transform: 



-| />00 -I /'OO 

1 [ s ) JO L K s ) JO 



(0|«£/9 



The modular transformation t = i-£f3 —> — - = i-Aig, or alternatively the Poisson summation formula, 
connects our Riemann Theta constant, which is a modular function of weight — g, to another Theta 
constant in the form: 



e 








I 



\mp) = -^=e 

1 l H> V4vr^ 








j AttI3> 
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The Mellin transform of this Poisson-inverted formula gives 

t / ^.-3 ^T^ -V 2 I3- 



(k (s) 



4vrr(s) Jo 
I 



/•OO _ _ 



4/3 



4vrr(s) 



T^ + 2 E 



,2s- 



n£Z/{0} 



n/ 



tfi ,(nJ) 



(16) 

(17) 



In formula (17), the Kelvin functions Ki_ (In) are transcendental entire functions of s (holomorphic for 
every s 6 C/oo, with an essential singularity at s = oo). Therefore, the poles (k (s) are the poles of the 
Euler r(s — g) function: 



113 5 



2j + l 



5 (-» ' r\ 1 c\ 3 5 



2' 2' 2' 2 



? 5 



JG^ 



The lengthy asymptotics, where the Hilbert spaces on the circle and the real line are traded, L 2 (S 1 ) 
L 2 (R), is easily read from formula (16). The Mellin transform is then applied to obtain: 



hK (P) 



1 

2^ 



cfte-W^) 



le~ v2 P 



V^P Cko(s) = ^ 



dk 



(k 2 + v 2 



^J r( s -J) 

/J^Jta-l T(s) 



(18) 



For later use it is also interesting to consider the high-temperature, small /3, asymptotics: 



e 









1±„2 



'4n/B' 



^ e ^ n ^0^o l + 0(e P) 



The idea of zeta function regularization now follows naturally: simply assign to the infinite vacuum 
energy a finite value provided by the spectral zeta function at a regular point in s £ C: 



A^)-fc,oH)-A^)) [s] .M^(^y CiCo( 



s) 



(19) 



where [i is a parameter of dimensions L 1 introduced to keep the dimensions of the regularized energy 
right. We stress that s = — g is a pole of this function. 

5.2 Kink spectral zeta and heat functions 

The spectral K-heat function with periodic boundary conditions at the finite length interval is: 



M/3) = TV L2(§1) e^ = 1 + J^ ^ + ^ T ' ''"' 



E 



(20) 



J=2 q„£Ker$(q,l,n) 

where the q n satisfy the transcendental spectral conditions 

1 2ir 

<5>(q,l,n)=q+j5(q)-—n = 

There is now a delicate point to be mentioned before of performing the Mellin transform. The zero mode 
oj 2 = in the spectrum of K introduces a spurious contribution in the spectral K-zeta function. In fact 



1 f'°° 1 

lim — - / dB I3 s ~ l e- £ P = lim — 
e^o F(s) J e^o e s 



(21) 
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valid only if Re(s) > 1, and it diverges in the e — > limit when it is defined. It is known, however, that 
the zero mode does not contribute to the one-loop quantum correction AE(cJ)k) [30] but only enters to 
second order in the h expansion. Fortunately, within the zeta function philosophy, the result (21) can 
be extended to the whole s-complex plane and it is if Re(s) < 0. Thus, we shall perform the Mellin 
transform of the i<f-heat function on the subspace orthogonal to the kernel of the K operator without 
any loss of information regarding our problem: 

1 Z" 00 
Tr i2(S1) e-^ =Tr i2{§1) e-^-l , ( K± ( s ) = d(3 /^Tr^^e"^ 

1 I s ) Jo 

b 1 1 

^( s ) = Er^+ E < a 2 +v 2 Y ■ ( 22 ) 

j=2 U 3 q n eKer<!>(q,l,n) Wn T ' 

At the I — > oo limit, the regularized K-heat and K-zeta, functions become: 

6-1 



h K ,ifi) = Tr LHm e-^ = Y^ e^' + e~^ + £ g (l + f£) e~^^ (23) 

, / n m /^±\~ s v^ 1 1 f°° dk f i d5\ 1 , ,, 

&x( S ) = IV I1(R) (^) =^ + +y ^, + _j^^^ . (24) 



J= 2 ~3 "b ^ 

Thus, the zeta function regularized kink Casimir energy reads as follows: 

„2 / „2\s+\ 
2^ ^ + ,2s, ,2s - TTTs + / O^St7fc2T^2V ' ' ^ 



A^(^)W = AE(<f> K )[s] - AE(^)[s] = ^ (j^J \ Ck± (s) - Ck (s)) 



1 4 2s f°° dkd5 i 

s7T2 
j=2 ~3 °b ^b 



The energy due to the one-loop mass counter-term 

AE2iM _ _M Wl)) £ ^*_ _ _M„ W) & J E 



/ie2, 



^n 2 + 1; 2 



is also regularized by the zeta function procedure: 

2 \ m dJ y s > 

Note that we regularize this quantity by means of Ck ( s + 1). This is necessary to compare residua at 
the same point of the s-complex plane. Finally, we obtain the zeta function regularized DHN formula: 

AE^k) = 

= M Urn i f4) S ' 2 f^{r^/ 3S " 1 (^4/ 3 )-^o(/3)) + ^lim^f^^„( 5 + l)) (26) 
2 s _> i \mij T(s) [J 2 i-K» I J 



W r ( s ) | ^ w ? s INfe S w2s V-oo 2vr 8k (k 2 + t; 2 )* + 2 w 2 ^ 1 
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5.2.1 The prototypes: sine-Gordon and X<p A kinks 
• In the sine-Gordon model we have: 



h Ko (P) 



I 



-P 



h K ±(P) = h Ko (P)-Erfc(y/p) 



<K (s) 



I T{s- 

In r( s ) 



(k^( s ) = (k (s) 



0r sT(s) 



We can now compute the limit of the regularized quantities that enter the one-loop correction 
formula to the kink mass: 



AE 



lWK) 



hx 

2m 2 
TiX 



lim 



M 



in- 



(C k ±(s)-Ck (s)) 



2m 2 



-l V m '' 

2 



s+i 



2 i ... I 
lim ' ' ' 



T{s + \] 



^JtTsT(s) 



In order to fall in the pole carefully, we look at the behaviour of this function in its neighborhood 
s 



-2 + e, and we take the e — > limit: 

nx 



AE^k) 



2y/Trm 2 

hX , 

o hm 

2irm 2 e^o 



lim ( — ' 

e->0 \m 

1 



Tie) 



J ra 



+ e) 



2 log -^--^(1) + ^) + 0(e) 
cm 2 



(27) 



where tp(z) 



Yr^Y is the digamma function. Simili modo, 



AE 2 



lK\ 



nx , /u\2s+ir(s + i 

— ^ hm — 

^ 1 \ TTJ / 



/?) 



/iA 



.1 Vm- 

2 



lim 



sm 2 s-*— - v '» 



r(s) Z-s-oo z 

M N2 S +ir( s + |) 



lim^Cxo(^ + l) + 0(^^: 



m 4 



r( s ) 



+ <w^S: 



;??, 



and 



AEo 



lK\ 



nx 



irm^ 



lim ( — 
e->o \m 



rfe) 



r(-i + 



+ ^ 2 S: 



m* 



o hm 



I + 2io g ii + ^(i)-^(-i) + 0( e ) 

£771 2 



+ 0(ft 2 ^ 



ttT 



The important point is that the renormalization of the zero-point energy performed by the sub- 
traction of Ck (— 5) stih leaves a divergence coming from the s = — ^ poles because the residua are 
different. The correction due to the mass renormalization counter-term also has a pole. The sum 
of these two contributions leaves a finite remainder and we end with the finite answer: 

AE(cj )K ) = AE 1 ( ( l )K ) + AE 2 ( ( t >K ) = -^ X 2 , (29) 

nm 

which is the same result as that obtained by cutting the number of modes in the DHN formula. In 
sum, the one-loop renormalized kink mass can be written in the form: 



£(V>l-loop) = %) + A%) 



8m 3 Km 



X 



7T 



+ o(h 2 



8m 3 



nr 



an m z mr ' 



8M 3 
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Note that if h^ > 8tt, the coupling constant is greater than the famous bound [34] where the 
quantum sine-Gordon model does not exist. We find that, beyond Coleman's bound, E(if)i_i oop ) < 
and the kink energy would be lower than the vacuum energy. The renormalized particle mass M is 
also pathological and becomes negative when M 3 < 0. 



In the \(f> model, we have 



h K i_ 



I 



2y/2ir 
I 



dke-^ k2+ V 



-4/3 



e -4/3 + e -3/3 Erf ^ _ Erfc 2 ^p 



Ck (s) 



J i r( 8 -|) 

/87F4S-5 r(a) 



Ck-l(«) = Ca"o( s ) + 



ttT(s) 



3 s+ 5 



p rl Q il 3 _ll 
1 2- r l[2' a ' 2' 2' 3J 



1 1 

4*s 



where 2-^1 [o, b, c, d] is the Gauss hypergeometric function. The power expansion of 2-P1, 



r(|) ^ (_i )J T(j + i)r( a + j + i) 
1^2^ 3*1-1 r(i + 1) 



3 l j\ 



tells us that, besides the poles of (k (s), Ck ± ( s ) has poles at s = — | + j, — % + j, — | + j, ■ ■ ■ 



nhn*+ l 2) j=0 

2 ~rJ' — 2 ~r^' — 2 

j G Z + U {0}; i.e., as in the sG soliton case, Ck {s) and ( K ±(s) share the same poles except s = \ 
but the residua in the A(^ 4 )2 model are increasingly different with increasing values of |Re s\. 

The kink Casimir energy is: 



AE 



1WK 



nx 

2m? 

hx 

2^2 



lim 



/' 



2 \ s 



.1 \ m* 

2 



(C K ^s) - Ck (s)) 



lim 

S^ 2 



~2 ) /=tV_n I 777T2^l[2' S + 2' 2' 



V^r( s ) V3 S+ 



13 ll 



1 1 

4^s 



and near the s = — I pole we find: 



AE 



1WK 



h\ In 

hm 



2 \ £ 



y^nm? e^-0 \m 2 J r( — 3+e 



r( £ ) 



A Ffl r 3 _ll J; 

Q£ 2- r H2' t ' 2' 3J / 1 . x .-I. 

6 (-i + e)4 2 



hx , 

~ hm 

2irm z £->0 



3 In 



//, 



777^ 



+ 2 + ln-- 2 ^[i,0,|,-i] + O( e ) 



where 2-^1 is the derivative of the Gauss hypergeometric function with respect to the second argu- 
ment. The regularized mass renormalization contribution to the kink energy is: 

Aa(*«) = -^i.m,(ii) 2 * +ir ^±iiii m %„(» + i) + o(tf4) 



7?) 



3hX , 

-== — hm 

/27T771 2 s->-± V "7 



r(s) i^-oo 1 



777,^ 



4 s +l r( s ) 



+ 0(fr 



A 2 



mr 



Near the s = — 5 pole we have: 

3frA 



AS 



2l<PX 



lill)( VVjr£(£L + 0(tf 



V2W ^0 V 777, 2 / r(-i+£ 



A?, 
4' 



/?? 



3/7A ,. 
k hm 

2-7T777 2 e^O 



-+ln^ ? -ln4 + ^(l)-^(-i) + C( £ ) 

e m z 2 



+ o{h 2 



mr 
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Again ,the sum of these two contributions leaves a finite remainder and we obtain: 

hm 3hm 



AE(i> K ) = AE^k) + AE 2 (<t> K ) 



2^6 vrv 7 ^ ' 



(30) 



the same answer as that discovered by regularizing the ultraviolet divergences with a cutoff in the 
number of modes. The one- loop renormalized kink mass is: 



£(Y>l-loop) = E{#k) + A£(^ K ) 



4 m A 



3^2A 



i-h- 



x 



rn- 



9 V3 

4^ ~ ~8~~ 



+ 0{ti 



2 A 2 



4 M 3 






3>/2A 








4tt 


1 
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Like in the sine-Gordon model we find a similar pathological situation: when h 

E(tpi-\ oop ) < and M 3 < 0. The real cubic root of this negative value of M 3 is negative such that 
M is also negative suggesting that this model is quantum mechanically unsound in this regime. 



6 The heat kernel expansion and the one-loop kink mass shift 

Both when we regularize the DHN formula by means of a cutoff in the number of modes and when we 
use the zeta function regularization procedure we must have a complete knowledge of the spectrum of 
the K Hessian operator governing the small kink fluctuations. For generic kinks this information is not 
available, and we need an alternative procedure to compute the one-loop mass shift other than the sG 
and \(j) A kink mass shifts. 

6.1 The asymptotic expansion of the kink heat function 

To perform this task, we shall rely on the asymptotic expansion of the heat-kernel expansion, which en- 
codes all the spectral information about the K differential operator. The kink spectral function r Yv L 2e~P K 
admits an integral kernel representation. If fk(x) are the eigenfunctions of K with eigenvalues w 2 (fc) la- 
beled by a continuous index k, the integral 

(31) 



K K (x,y;P) = jdkf* k (y)f k (x)e-^W 
is the K-heat equation kernel. The K-heat function is thus: 

/oo 
dxK K (x,x;P) 
-oo 

The .ftTo-heat equation kernel is found immediately: 

K Ko (x,y;fi) 



(32) 



V^P 



Rl ,2 _ (*-«) 



because in this case: fk(x) = ——e and w (k) = k + v , k £ R. 



To estimate Kx(x,y;f3) in the cases where the eigenfunctions of K are not known we shall use the 
fact that this integral kernel is the fundamental solution of the K-heat equation: the solution of 



d d 2 



+ v 2 + V{x) 



d(3 dx 2 
with a (5-source on the diagonal of M 2 at infinite temperature: 

K K (x,y;0) = 6(x-y) . 



K K (x,y;[3) = 



(33) 



(34) 
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We shall assume the convolution ansatz 

K K (x,y;f3) = K Ko (x,y;(3)C(x,y;f3) (35) 

in order to benefit from the knowledge of KK (x,y; /3). Plugging (35) into the if- heat equation (33) we 
find the transfer equation for C(x, y; /?): 

h +x -^i-^ +v(x) ) C(x - v ^ )=a ■ <36) 

Moreover, the infinite temperature condition (34) prompts the condition 

C(x,y;0) = l . (37) 

The idea is now to solve (36) by means of a series expansion. Firstly, C(x, y; /3) is expanded around the 
point (3 = as a power series in f3: 

oo 

C(x,y;(3) = Y,c n (x,y)P n . (38) 

n=0 

The zero-order density is fixed by the the infinite temperature condition (37): 

c (x,y) = l . (39) 

Secondly, (38) is plugged into (36), which trades this PDE for the recurrence relations between the 
densities c n (x,y) and their first- and second-order derivatives: 

(n + 1) c w+1 (x, y) + (x - y) dCn+ ^ V) + V(x)cn(x, y) = ^f^ v) . (40) 

The solution of (40) provides a series expansion of the integral kernel in the form: Kx(x,y; (3) = 
Kk (x, y; j3) J2^=o c n(x, y) j3 n . In turn, the if -heat function is also expanded in powers of (3 

e -/3v 2 °° 



h K ((3) = - / ==Y, c ^ K ^ n . ( 41 ) 

n=0 

where the c n {K) are the Seeley coefficients of the if- heat function high-temperature expansion: 

/oo 
dxc n (x,x) . (42) 

-oo 

The formula (42) above involves the limit of the densities c n (x, y) as y approaches x, i.e., 

c n (x,x) = limc n (a;,y) . 

Moreover, derivatives of the densities enter the recurrence relations (40). It is thus convenient to define 
the y — $■ x limit of the derivatives 

^C n (x) = hm ^%^ (43) 

y^>x OX K 

such that 

c n (x,x) = ^C n (x) . 
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The infinite temperature conditions (39) determine the densities ' 'Co(x): 

^C (x) = lim^ = S k0 

Taking the fc-th derivative of expression (40) with respect to x and later taking the limit y — > x, we 
obtain recurrence relations between all the densities ' >C n (x): 



(k) 



C n (x) 



(44) 



affording us the identification of ' >C n (x) in a recursive way from ^ >Cq(x) = 5 . Note that the compu- 
tation of the densities ^ k >C n {x) requires the densities (°'C n -i(x), ^ 1 'C n -i(x), . . ., ^ k+2 'C n -i(x), which in 
turn are determined from the densities (°)C n _ 2 (x), ( 1 >C n -2(x), . . . ' fc+4 )C n _ 2 (x) and so on until we reach 
the densities ^Cq(x) = 5 k0 . Therefore, the evaluation of the density ( k 'C n (x) is based on the densities 
that appear in the following pyramid: 

(0) C n (x) 

(°)Cn_i(x) (1) C n _!(x) ^C n ^(x) 

(°)C„_ 2 (x) WC n . 2 (x) ( 2 )C n _ 2 (x) ( 3 )C n _ 2 (x) ( 4 )C n _ 2 (x) 

^Co(x) (1) C (x) {2) C {x) (3) C (x) {A) C (x) ... ^ 2n -^C G {x) (2n) C (x) , 

i.e., (n + l) 2 densities ^>C m {x) must be evaluated previously. 

We list the first four densities c n (x, x) for the kink fluctuation operator in a general (1+1) dimensional 
one-component scalar field theoretical model: 

co(x,x) = ®C (x) = 1 

ci(x,x) = (°)Ci(x) = -V{x) 

c 2 (x,x) = Mc 2 (x) = ~\^ + \{V{x)f 

, ^ m„ , n 1 d 4 V 1 , ,d 2 V 1 dVdV \. ,3 

c 3 (x,x) = (0) C 3 (x) = -rr^- + -V(x)—^r H — - V(xY . 

y ' y ' 60 <9x 4 6 v ' dx 2 12 dx dx 6 y ' 

There is an interesting point about these densities: they are the infinite conserved charges of the KdV 

equation, see e.g. [3(3] . Consider the family of differential operators 

K(/3)=-^+V 2 +V(x,P), 

where the family of potentials solve the KdV equation: 

dV_ dV d 3 V _ 

8(3 dx dx 3 

The /3-evolution of K{{3) can be expressed in the Lax pair form 

§ + [*M1-0 , MW ) = 4^ + 3 (vl + lv) + B im 



such that it is iso-spectral: 

K((3)=g(f3)K(0)g- 1 (P), M{P) = %~ x {fi) . 
The densities codify the spectrum of K. Ergo, 

dx e -j= V c n (x, x- f3)(3 n => c n (x, x- (3) jL /(/3) 
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6.2 The asymptotic high-temperature expansion of the quantum correction 

The K-heat function h K ±((3), basic in the computation of the one-loop mass shift via zeta function 
regularization, admits the series expansion: 



lim 

l— >oo 



-fiv 2 -j3v 2 -fiv 2 °° 

h K ±[fi) = e —=^{K)ii-\ + C1 (K)^ + e -= Y,c n {K)f3 n -l - 1 

V47T V47T V47T *- 1 

n=2 

The first two Seeley coefficients are: 

/"* 

cq{K) = I co(x,x)dx = l and c\{K) 

The zero-point energy renormalization 

h K x(0) - h Ko (P) = ^La(K)ph + e —= ^ Cn (^)/3"-2 - i 

^ 47r n=2 



(45) 



ci(x,a;)(ix 



F(z)dx = - (V(x)) 



-/3v 2 



-in 



(46) 



cancels the (divergent at the I = oo limit) contribution to the heat function of the cq(K) coefficient. 

Another crucial cancelation, known as the heat kernel renormalization criterion, see e.g. [16], occurs 
between the contribution of the c\(K) to the one-loop mass shift and the mass renormalization energy: 
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pole are canceled and we finish with the following one-loop 



Therefore, the divergences in the s — 

mass shift formula derived from the high-temperature expansion of the K- and K^-heat functions: 
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A*) = f lim , 

2 s ^_i r(s) j 



d/3/3 



8-1 



-I3v 2 



Y,Cn{K)p n -^-l 



n=2 



(47) 



7 A computational treatment of the one-loop kink mass shift formula 

7.1 The computational formula 

In this Section we shall convert the one-loop kink mass shift formula (47) into a formula that will afford 
us a precise estimation of the quantum correction in a general model simply by using a desk computer. 
In this computational scenario the series (47) is first truncated at a certain order N. Before doing all this 
we note the following behaviour at high and low temperatures of the difference of the K- and KQ-heat 
functions: 



]im(h K x(P)-hK (J3)) 

/S-K> 



-1 



lim (h K x(P)-h Ko (P))=0 ■ 



We shall denote as Sk(/3, N) the truncation of the series expansion (46) of this quantity: 



Sk((3,N) 



-/3v 2 



-pv 



2 N 



4tt 



Cl (inM + ^=]Tc n (K)/r-i-i . 



(4* 



n=2 
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This truncation required for computational effectiveness involves a troublesome effect. The asymptotic 
behaviours of the series (46) and the truncated series (48) for Hk{P) — hK (P) do not coincide at low 
temperature. It is easy to check that: 



\jmS K (P,N) 



lim S K (P,N) = -1 

/3— >oo 



As expected, the behaviour of h K ± {(3) — hx (/?) and Sk(P, N) fits extremely well in the high-temperature 
regime, small values of /3. If N is large, the identification of the series and the partial sum can be continued 
up to a value /3 = Po(N) such that h K ±(f3o(N)) — hK {Po(N)) is close to zero. The convergence of the 
partial sums to the series implies that limTv-s-oo Po(N) = oo. This is depicted in Figure 2 for the A</> 4 
model and several values of N. Replacing h K ±(/3) — h,K (P) by the partial sum Sk((3, N) in formula (47) 





Figure 2: A</> 4 kink: function h K ±(f3) — fiK {/3) (left) partial sums Sk{(3, N) for growing values of N from 
the bottom (red) to the top (violet) line (right), both depicted as functions of (3. 

we obtain: 
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2 s ^_i r(s) j 



d(3 



-fiv 



2 N 



'Air 



X>( x )£ 



s+n- 
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s-1 
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+ R(Po,N) 



where R(/3q,N) is a negligible contribution such that: 

lim R(Pq(N),N) = . 

N— >oo 

The next steps are to perform the /3 integration and to take the s — > — \ limit: 



AE^k-N) 



M lim i 



-§r( s ) 



1 N , 

= J2cn(K)(v 2 )*- n - s 1 [-- + n + s,Pov 2 



ft 



n=1 



2\/4^ 



1 N 

Y J c n(K)(v 2 ) l - n 7[n-l,(3 v 2 } + 

47T — ' 



n=2 



% 



where j[— \ + n + s,Pqv 2 ] are incomplete Euler Gamma functions. Finally, 



AE^k; N) = h 7 2 d 



1 N 
--Y^CniK^f-^in-l^ov 2 



n=2 



2V^ 



(49) 



is the basic formula for our computational approach. 



7.2 Constructing a symbolic algorithm 



We are now ready to construct a program by means of a symbolic software application such as Mathe- 
matica. The key process is to implement formula (49) in this algorithm. In order to obtain an estimation 
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of AE((f>K)N) we must compute the c n (K), n = 1, 2, . . . , N, Seeley coefficients. This computation re- 
quires the integration of the c n (x,x), n = 1,2, ...,N, densities over the whole real line. Previously, 
we need to identify (JV + l) 2 densities ( >C n (x) using the recurrence relations (44) in order to evaluate 
cn(x,x) = (°)CV(x). For instance, the calculation of AE(c/)k;20) demands the computation of 441 den- 
sities of this type. Once the Seeley coefficients have been identified we apply formula (49) in a direct 
way. All these tasks are divided into three subroutines described below: 

• Calculation of the c n (x,x) densities. 

A key point in our symbolic program is the possibility of expressing all the densities as functions of 
the 4>k field without the need to use their explicit dependence on x. Once we identify the ^ k 'C n -i(x) 
as functions of 4>k(x), the next-order densities ( k 'C n (x) can also be written as functions of the field. 
The only elements in the recurrence relations (44) that we need to determine in terms of (ftx are 
the derivatives -^ . To do this we recall that V(4>k) = ^ti[4>k] — v 2 and we use the first order 



ODE (7): 



dV dV d<j> K d 3 U r—— 

Ix- = ^ K] ^ = W [4>KW [M 

d 2 V fd 4 U 



dx 2 



[faWwifi 



K 




d 4 U dU 3 

= V^bM+fi 

Higher than second-order derivatives -j-f,j > 3, can be expressed as functions of (J)k by following 
the same pattern. A symbolic program easily manages this computational strategy in each step 
from n — 1 to n in order to calculate the ^ k 'C n (4>K) densities in an optimum way. For instance, we 
find that: 

d 2 U 

Cl(x,x) = V 2 - -Qlpl^Kix)} 

t ^ 1 ( 2 d2U ^ , mV ldU ^ < M d3U ^ < M 1 TT^ , M d * U ^ < M 

C2{X,X) = - \V ~ -QJp[<t>K{x)]j -Q^[4>K{x)]^[4> K (x)]--U[<t>K{x)]^[<l>K{x)] 

and so on. The following Mathematica code 

densitycoef f icients [potential., vacuuml_, vacmun2_, nmax_] := Module [{var4, 
var5, var6, tomax, dl, v, vO, oper, f6, x7, coef icientes , alfa, coa, k8, co} , 
(var4[phl_] = potential/ . {y -> phi}; var6 [phl_] = Simplify [PowerExpand [Sqrt [2 
var4[phl] ] ] ] ; var5[phl_] = Sign [var6 [1/2 (vacuuml+vacuum2) ] ] var6 [phi] ; coeficientes 
= {}; v[x_] = Simplify [(D[var4 [phi] ,{phl,2}]) /.{phi -> phi [x] }] ; vO[x_] = 
Simplify [(D[var4 [phi] ,{phl ,2}] )/ . {phi -> vacuuml}] ; dl [f un_] := Simplify [(D [fun, x] ) 
/. { phl'[x]-> var5[phl[x]]}] ; oper[fu8_, nlj := Simplify [Nest [f 6, x7, nl] /. {f6 
-> dl , x7 -> fu8}] ; tomax = 2 nmax; For[alfa = 0, alfa < tomax + 0.5, coa[0, alfa] 
= 0; alf a++] ; coa[0, 0] = 1; co[k5_, alfaj := Simplify [l/(k5 + alfa) (coa[k5 - 1, 
alfa + 2] - ]Cr5= a Binomial[ alfa, r5] oper[v[x] - vO [x] , r5] coa[k5 - 1, alfa - r5] 
) ]; For [k8 = 1 , k8 < nmax + 0.5, tomax = tomax - 2; For [alfa = 0, alfa < tomax + 
0.5, coa[k8, alfa] = co[k8, alfa]; If [alf a == 0, coeficientes = Append [coef icientes, 
coa[k8, 0]]]; alf a++] ; k8++] ; Return [coef icientes] )] ; 

defines the module densitycoef f icients [potential., vacuuml., vacuum2_, nmax _] that is ca- 
pable of performing this work. The arguments of this computational function are potential, the 
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U(y) potential written by prescription as a function of the y variable, vacuuml and vacuum2, the 
two vacua connected by the kink solution in increasing order, and nmax, the iV-order truncation 
chosen in the computation of AE((/)k', N) . 

Calculation of the Seeley coefficients. 

We proceed to the calculation of the c n (K) Seeley coefficients by integrating the c n (x,x) densities 
over the whole real line. A novel strategy to perform such integrations will reduce the main 
difficulty; namely, to cope symbolically or numerically with infinite integration domains. Recall 
that the densities c n (x,x) are functions of the kink field 4>k'- c n (x,x) = / n (0x)- The change of 
variables z = 4>k{x) is a bijective function from [—00,00] to \z\ = <ft( l > , 22 = <ft ] such that: 

f°° f Z2 fn(z) 

c n (K) = / f n ((j) K (x))dx = / — = dz 

J-oa Jzi \/2U{z) 

because dz = —^-dx = y/2U(<pK)dx. This manoeuvre avoids the problem of evaluating an 
integration over an infinite domain and explicitly shows the topological character of the See- 
ley coefficients, which only depend on the values of certain functions of the field -let us call it 
generalized superpotentials- at the vacuum points. The following Mathematica module called 
seeleycoeff icients [potential., vacuuml., vacuum2_, nmax_] 

seeleycoef f icients [potential., vacuuml_, vacuum2_, nmax_] := Module [{coef , 
f, fl, f2, a = {}, k8} , (coef = densitycoeff icients [potential, vacuuml, 
vacuum2, nmax]; fl[y_] = Simplify [PowerExpand [ Sqrt [2 potential] ]] ; f2[y_] 
= Sign[fl[l/2 (vacuuml+vacuum2)]] fl[y]; For [k8 = 1 , k8 < nmax + 0.5, f[y_] = 
Simplify [(coef [ [x8] ] / . {ph[x] ->y} )/f2[y]]; a = Append [a, Integrate [f [y] , {y, 
vacuuml, vacuum2}]]; k8++] ; Return[a])]; 

does the trick. Here again the arguments of this function are the potential U(y), the vacuum points 
<jft-' and 4>( 2 > and the truncation order JV. Note one remarkable fact: there is no need to know the 
explicit expression of the kink solution !! The potential U(cp) and the vacua connected by the kink 
encode all the necessary information in this computation. 

From a computational point of view this subroutine is vital for the precision of the final result. 
Some dysfunctional behavior may arise from this algorithm when the value of the truncation order 
N is large. This may happen when the densities c n (x, x) include such large/small factors than some 
loss of precision occurs. This effect can be detected by checking the sequence of Seeley coefficients. 
Theoretically these coefficients tend to vanish when the index n increases. Thus, if there is some 
anomaly in the Seeley coefficient sequence starting in a given no the quantum correction estimated 
as an asymptotic series can be trusted only up to a truncation order N = uq. 

Estimation of the quantum correction. 

It remains to optimize the value of the /3o(iV) choice, related to the asymptotic behaviour of the 
partial sums Sk(/3,N). Firstly, we find a point (3 at which Sk(P,N) becomes approximately —1 
and, second, we move back until we obtain a value /3o such that Sk(P,N) almost vanishes. With 
this choice of optimum j3o we obtain the quantum correction by applying (49), and the task is 
accomplished by the subroutine quantumcorrection [potentials vacuuml_, vacuum2_, nmax_] 
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quantumcorrection[potential_, vacuuml_, vacuum2_, nmaxj := Module [{var4, 
vhessl, heat, a, bO, corr}, (var4[phl_] = potential /. {y -> phi}; vhessl = 
Simplify [(D[var4 [phi] , {phi ,2}] )/. {phi -> vacuuml}] ; a = seeleycoef f icients [potential, 
vacuuml, vacuum2, nmax] ; heat [b_] = E~{-b vhessl}/(2 Sqrt [Pi] ) (J2n=T (a[[n]] b n_ 2) - 
1; bO = 10; While [heat [bO] > -0.9, bO = bO + 1] ; For [h = . 1 , h > 0.001, h = h/2, bO = 
bO + 2 h; While [Abs [heat [bO]] >= Abs [heat [bO - h] ] , bO = bO - h]]; corr = -1/(2 Sqrt [Pi 
bO]) - 1/(8 Pi) (En™ (a[[n]] vhessl~ n+1 (Gamma[n - 1] - Gamma[n - 1, bO vhessl]) )); 
Return [corr] )] ; 

The above Mathematica code characterizes the application that automates the computation of the quan- 
tum correction to the kink mass in one-component scalar field theory. 

The KinkMassQuantumCorrectionJLite.nb file containing this Mathematica code can be download at 
the web page http://campus.usal.es/~mpg/General/Mathematicatools. We recommend this option in 
order to avoid transcription errors in the code. 

7.3 Testing the algorithm in the sine-Gordon and A</> 4 models 

7.3.1 The algorithm and the sine-Gordon kink 

The exact one-loop sG-kink mass shift, calculated by means of the DHN formula regularized either by 
a cutoff or by the zeta function procedure, is: AE^q \4>k) = —hm(-) ~ — 0.31831/im. The result for 
this magnitude provided by the asymptotic approach is accessible in the symbolic algorithm. Recall that 
r )\ = vn and that the dimensionless potential is U((f>) = 1 — cos</>. By running the KinkMassQuantumCor- 
rection_Lite.nb file and then executing the command quantumcorrection[l — cos (y),0,27r, nmax] with 
nmax = 10, 20, 30 the following results are obtained: 

quantumcorrection[l — cos(y) ,0 ,2ir , 10] =-0.3191 13 
quantumcorrection [1 — cos(y) ,0 ,2ir , 20] =-0 . 318321 
quantumcorrection [1 — cos(y) ,0 ,2ir , 30] =-0 .318310 

The relative errors are: 1) 0.252263%, truncation at 10 terms, 2) 0.00345755%, truncation at 20 terms, 
and 3) 0.0000600568%, truncation at 30 terms. The final result fits in extraordinarily well with the exact 
one-loop sG-kink mass shift and it is closer for higher truncation order. 

The symbolic algorithm is prepared such that it will provide not only the final outcome but also offer 
the intermediate quantities needed in the computation. For instance, the sG-kink Seeley coefficients are 
accessible by executing the command seeleycoef f icients [1 — cos(y) ,-1,1, 30] . See the list of the 30 
first coefficients {c n (K)}^ =1 : 

{4., 2.66667, 1.06667, 0.304762, 0.0677249, 0.0123136, 0.0018944, 0.000252587, 0.0000297161, 
3.12801 • 10- 6 , 2.97906 • 10~ 7 , 2.59049 • 10~ 8 , 2.07239 • 10~ 9 , 1.5351 • 10~ 10 , 1.05869 • 10" 11 , 
6.83027 • 10~ 13 , 4.13956 • 10~ 14 , 2.36546 • 10" 15 , 1.27863 • 10" 16 , 6.55706 • 10" 18 , 3.19857 • 10" 19 , . 
1.48771 • 10" 20 , 6.61203 • 10" 22 , 2.81363 • 10~ 23 , 1.14842 • 10" 24 , 4.50361 • 10" 26 , 1.69947 • 10" 27 , 
6.17991 • 10~ 29 , 2.16839 • 10~ 30 , 7.35047 • 10~ 32 } 

The command densitycoeff icients [1 — cos(y) ,-1, 1,N] unveils the densities c n (x,x) needed in the 
definition of the Seeley coefficients (42). The general form of these densities is: 

f \ 2 " ■ 2 <t>K 2"- 1 

c " (x ' x) = (2^iyyy sin t = (2^i)n c/( ^ ) 

and, henceforth, 

. . 2 n f 27T zdz 2 1+n 

c n( K ) = l^r T\77 / sm 



(2n-l)!! Jo 2 2 (2n-l)!! 
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7.3.2 The algorithm and the A0 4 kink 

The exact DHN one-loop A0 4 -kink mass shift is: 



AK, 



km 



'>K 



y/2 V2\/3 vr 



-0.666255 



Km 

V2 



Now 7J = -7= and the dimensionless potential is U(cf)) = ^ (</> 2 — I) 2 - Therefore the expected dimensionless 
quantum correction to the kink is —0.666255. Again, running the KinkMassQuantumCorrectionJLite.nb 
file on a Mathematica platform and then executing later the command quantumcorrectiont^y 2 — 
l) 2 , -1,1, 10] we obtain the estimation of the one-loop quantum correction to the kink mass in the 
asymptotic approach truncated to 10 terms. The value of N must be chosen adequately so that N 
will be large enough to obtain a precise answer but not so large that the process becomes very much 
time consuming. The processing time is very sensitive to the complexity of the potential U(y); note 
that the symbolic software algorithm must manipulate and simplify the function U(y) and its successive 
derivatives. Trying three increasing values of N, we obtain: 

quantumcorrectiont^y 2 - 1) 2 ,-1, 1, 10] =-0.665894 
quantumcorrectiont^y 2 - 1) 2 ,-1, 1,20] =-0.666241 
quantumcorrectiont^y 2 - 1) 2 ,-1, 1,30] =-0.666254 

The relative errors are: 1) 0.054112%, truncation at 10 terms, 2) 0.00202985% truncation at 20 terms, 
and 3) 0.0000786421%, truncation at 30 terms. These results guarantee the precision of the method. The 
command seeleycoeff icients [^(y 2 — l) 2 , -1,1,30] prompts the list of the 30 first Seeley coefficients 
c n (K), n = 1, . . . , 30, for the A</> 4 -kink. 

{12., 24., 35.2, 39.3143, 34.7429, 25.2306, 15.5208, 8.27702, 3.89498, 1.63998, 0.624754, 0.217306, 
0.0695378, 0.0206038, 0.00568381, 0.00146679, 0.000355585, 0.0000812766, 0.0000175733, 
3.60478 • 10" 6 , 7.03372 • 10~ 7 , 1.3086 • 10~ 7 , 2.3264 • 10~ 8 , 3.95983 • 10~ 9 , 6.46503 • 10" 10 , 
1.01412 • 10~ 10 , 1.53075 • 10" 11 , 2.22655 • 10~ 12 , 3.12498 • 10~ 13 , 4.23726 • 10~ 14 } 

The command densitycoeff icients [^(y 2 — 1) 2 ,-1, 1,N] in turn provides us with the densities c n (x,x): 
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leading to the general form of the A(/> 4 -kink Seeley coefficients: 
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Figure 3: Seeley coefficients depicted as bar charts for the <fi model (left) and the sine-Gordon model (right). The 
one-loop shifts are also shown in red 

We now emphasize the strong points in these new calculations. 
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• In comparison with the primitive symbolic algorithm used in Reference [28], the present Mathemat- 
ica program has been extraordinarily refined. The error committed with respect to the exact result 
in the one-loop sG-kink mass shift truncated at N = 10 was 6.00 percent and the relative error 
committed in the present calculation is 0.25 percent. In the X(p 4 kink case, the past and present 
errors are 0.07 and 0.05 percent. It is the better fitting of the parameter /3q in the sG-kink that 
helps to greatly improve the approximation. Also, integrating in the field space rather than in x is 
more beneficial in the sG kink because the exponential tails of the kink solution decay more slowly 
to the vacuum values than their counterparts in the Xcp 4 kink. 

• Skipping the integration over infinite domains also helps to compute many more coefficients. A 
frequent criticism of the Gilkey-de Witt procedure, the difficulty of computing many coefficients, is 
thus circumvented in this relatively easy problem. Needless to say, the errors are hugely reduced. 

• In the basic sG and A(/> 4 kinks, the functional dependence on the kink profile of the Seeley densities 
is independent of n. This is not so, as we shall see, in other cases, hindering the calculation of the 
Seeley coefficients. 

8 Other models: one-loop corrections to the mass of exotic kinks 

It is possible, in principle, to apply the algorithm to compute the one-loop kink mass shifts in any other 
one-component scalar field theory model. The inputs that must be introduced in the algorithm are very 
few: the (dimensionless) potential supporting kinks; the nearest neighbors vacuum points connected by 
the kink that must belong to a single point in the vacuum moduli space, and the truncation order. We 
stress that the choice of an excessively large value of N can finally cause a loss of precision, a feature 
characteristic of asymptotic series. 

8.1 One-loop kink mass corrections in distinguished models 

In Table 1 we assemble the inputs necessary to determine the kink masses up to one-loop order as well 
as the one-loop kink mass shifts and truncation order TV" in several models discussed in the Literature. 
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Table 1: The one- loop kink mass shift in several one-component scalar field theory models. 

The first model is simply the sine-Gordon model with a re-scaling of the scalar field: (f) —> 2(f). In 
this re-scaled sine-Gordon model the distance between two nearest neighbors vacuum points is only n. 
Therefore, we have a very good test to the algorithm because a simple change of variables unveils the exact 
relationship between the re-scaled sG kink and the proper sG kink mass shifts. One easily concludes that 
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the one-loop mass shift of the re-scaled sG kink must be twice (in absolute value) the one- loop correction 
of the sG kink mass because the Hessian operator on the re-scaled sG kink is four times the Hessian 
operator on the sG kink. The estimation shown in the Table is in perfect concordance with this fact. 
Note, however, that the classical mass behaves in the opposite way: the classical mass of the re-scaled 
sG kink is half the classical mass of the proper sG kink. 

The following entries in Table 1 come from other models. The unifying feature of all these other 
models is the impossibility of applying the DHN formula because of the lack of information about the 
spectra of the Hessians. Therefore, all the one- loop corrections given in the Table above are computed 
by means of the asymptotic formula. The potential in the second row is the famous Razavy potential 
of molecular physics, see [39]. Since the potential of the physical pendulum gives birth to the (1+1)- 
dimensional sine-Gordon field theory, we use the Razavy potential of molecular physics to build another 
(l+l)-dimensional field theoretical model. The model in the third row is akin to the </> 6 Lohe/Khare 
model [41], [42]. The difference is that in this case the minimum at (j) = becomes a maximum and the 
vacuum orbit has only two points, see Table 1. In rows 4 and 5 we respectively show pure cj) S and </> 16 
models. Note that in these cases only implicit expressions of the kink solutions are available. Thus, there 
are no analytic expressions of the Hessians on these kinks available. The ability to compute the one-loop 
mass shift in this situation is a very strong aspect of the asymptotic approach: one is able to compute 
spectral information even if the analytical expression of the differential operator is not known !!. The 
required input is only the potential function U((j>). The model in the last row is very special. It has been 
derived from the deformation method, see [40], and its moduli vacuum space has two points. There are 
thus six kink/antikink solutions 



/arccos(tanhx) + (j 

cos 
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j = l,2,-..,6 



but only the one shown in Table 2 (and its antikink) interpolates between vacua belonging to the same 
point of the vacuum moduli space. In Table 2 we offer the pertinent data of all these models. 
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Table 2: Kink solutions, superpotentials, and particle masses v 2 of the one-component scalar field theory 
models chosen above. 

We skip giving the specific parameters m^ and jd in each model, but its identification is not difficult 
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in comparison with the physical parameters chosen in the Literature. 

8.2 One-loop kink mass corrections in two one-parametric families of models 

We finally consider two families of real scalar field theory models with trigonometric and polynomial 
potentials depending on a real parameter a. The family of potentials in the first case is a family of double 
sine-Gordon potentials where the potential U^q is continuously deformed to the potential J7q . This 
example is a new stress test for our method because the quantum correction must interpolate between 
two known values. The second family of models is formed by a one-parametric family of generalized 
6 models with only two vacuum points that can be considered as a deformation of the </> 6 Lohe/Khare 
model, where a third vacuum arises, see [41]-[42]. 

8.2.1 One- loop kink mass shifts in a family of double sine-Gordon models 

The first family of models is characterized by the following family of potentials in (1): 

U((f);a) = 1 - (1 - a)coscj)- acos(2<f>) , 0<o<l,a£l . 

For fixed a, such that < a < 1, it is the celebrated double sine-Gordon model that has been studied by 
many authors. We specifically mention Reference [37] as the starting point for our calculations of one-loop 
kink mass shifts in the double sine-Gordon models; closely related models are considered in Reference 
[38]. Given a, £/(</>; a) is a non-negative function that, when a varies in the above range, interpolates 
between the sine-Gordon, a = 0, and the re-scaled sine-Gordon, a = 1, models: 
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For all the parameter values such that < a < 1 the vacuum orbit is Ai = {2nk}, fceZ. When a 
becomes greater than but is strictly less than 1, new relative minima (false vacua) emerge between each 
consecutive pair of vacua. At a = 1 the new minima become absolute minima, the vacuum points of the 
re-scaled sine-Gordon model, see Figure 4 (left), such that the vacuum orbit becomes M. = {irk}, k £ Z. 
The kink solution of each member of the family with a strictly less than 1 is: 
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yi — asinh(\/l + Sax) yl + 3a 

see again Figure 4 (middle left). The classical energy can be calculated from the superpotential 
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We can check that when a = the kink solution is 4>k(x',0) = 4 arctan e x : the kink of the sine-Gordon 
model. Its density energy is localized around one point. As the parameter a approaches a value of 1 the 
solution (50) is formed by two separate kinks somehow related to the kink of the re-scaled sine-Gordon 
model. In this situation, the energy density shows two lumps of energy, see Figure 4 (middle right). 
The splitting of the two lumps is only appreciable near the value a = 1. Note that for a = 0.9999 the 
distance between the two lumps is less than eight spatial units, see Figure 4 (middle left). This is also 
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manifest in the classical energy, see Figure 4 (right): if a = the classical energy is 8 as it should be. 
For a = 1 — e we find almost 8 again, meaning that we meet two re-scaled sG kinks well apart from each 
other. The consequence is that we must expect very close values to a = 1 in order to observe the effect 
of the splitting in the kink mass quantum corrections. For completeness we give the particle masses and 





Figure 4: Graphical representation of the family of potentials fora = 0, 0.25, 0.5, 0.75, 1 (left), kink solution profiles 
(middle left), kink energy densities (middle right) for a — 0, 0.25, 0.5, 0.75, 0.99, 0.9999 in several double sine-Gordon 
models. The kink classical energy as a function of a (right). 



the quantum potential wells: 
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which are not, in the new framework of this paper, necessary to compute the one- loop mass shifts. We 
shall work the asymptotic approach to compute the one-loop quantum correction to the masses of the 
kinks (50) for several values of the parameter a. The Seeley coefficients are depicted for the values 
a = 0.1, 0.2, 0.3, . . . , 0.9 in Figure 5. 




Figure 5: Seeley coefficients in the family of double sine-Gordon models for a — 0.1, 0.2, 0.3, . . . , 0.9. The lines 
join the values of the coefficients associated with a specific value of a. The coefficients increase with increasing a. 

The quantum corrections to the kink masses are given in Table 3, where the truncation order is also 
depicted. All these data are graphically represented in the attached Figure 5. Bearing in mind that at 
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Table 3: The value of the quantum corrections to the kink masses in the family of double sine-Gordon 
models for several values of the parameter a (left) and its graphical representation (right). 
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the limit a — > 1 solution (50) is formed by two solitons of the re-scaled sine-Gordon model, it is tempting 
to interpret the quantum correction AE at this limit as the double of the quantum correction of the 
soliton in the model mentioned. In particular, we recall that A£" sG = — 0.636894/im, see the results in 
Table 1, and that the classical mass is twice the classical mass of the re-scaled sine-Gordon kink. We 
observe an approximated behaviour like this in Table 3, although the accumulation near the point a = 1 
does not allow us to reproduce the answer with complete precision. 

8.2.2 One-loop kink mass shifts in a family of generalized 6 models 

The family of potentials in the action (1) chosen now is 

l^;a) = ^ 2 + a 2 )(<A 2 -l) 2 
The vacuum orbit is Ai = {—1,1} and the kink solution that connects these vacuum points are 
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From the family of superpotentials 
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we read the classical energy of each member of the kink family: 

E{cp K {a)) 



i (( ,_ 2)V ^77 + ^ + 4K ( log ^±l) 



In Figure 6 some members of the family of potentials, some kink solutions, and their energy densities for 
several values of the parameter a are depicted. A plot of the classical energy as a function of a is also 
offered. The particle masses and the Schrodinger potential wells are in this case 





Figure 6: Graphical representation of potentials (left), kink solution profiles (middle left), kink energy densities 
(middle right) for a = 0.001, 0.01, 0.25, 0.5, 0.75, 1, 1.25, 1.5 in some generalized 4> e models. The kink classical energy 
as a function of a (right). 
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Now, the computational approach unveils the one-loop kink mass shifts shown in Table 4 for several 
values of the parameter a, with a truncation of N = 11 terms. All these data are graphically represented 
in the attached Figure in Table 4. 
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Table 4: Quantum corrections to the kink masses in some members of the generalized 
several values of the parameter a (left) and its graphical representation (right). 



models for 



We remark an interesting behavior. At the a — > limit the potential reduces to U(<p, 0) = ^^i^? — I) 2 , 
i.e., the pure (j) G Lohe/Khare model; a new vacuum point <fi = arises, and the vacuum orbit becomes: 
M. = { — 1,0, 1}. The lim a ->o cj>k(x; a) goes to a configuration where the 4>k{o) kink solution splits into 
two kinks, one connecting the vacua <j> = — 1 and <f> = and the other connecting <p = and = 1, each 
one close to the kink of the Lohe/Khare model, see Figure 6. Note that as the parameter a vanishes 
solutions (51) involve two separated lumps of energy. We observe in Table 4 that the quantum correction 
to these particular kinks tends to infinity. The reason for this stems from the fact that the vacua that 



are connected by this kink give rise to different meson masses: 
References therein to see a treatment of this problem. 
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9 Summary and outlook 



In this paper and several earlier ones, [28], [33], [6], [7], we have elaborated a computational procedure to 
calculate one-loop kink mass shifts summarized in formula (49) supplemented with a very sophisticated 
use of symbolic Mathematica software. The derivation of this formula is based on the high-temperature 
heat-kernel expansion plus the Mellin transform, thus ending in some series of asymptotic type. The 
old DHN approach to this problem led, after zero point renormalization by subtraction of the vacuum 
from the kink energy mode-by-mode, to formula (12) which can be regularized either by a cutoff in the 
number of modes (13) or by means of the generalized zeta function procedure: (26). Except for the 
sine-Gordon and \<j) kinks these formulas are useless from a computational point of view because of 
the lack of information about the spectra of the differential operators governing the fluctuations around 
other interesting kinks. More precisely, what we have achieved in this work is some conceptual advance 
in the understanding of formula (49) derived (and exploited) in the papers quoted above that hugely 
augmented the precision of the symbolic algorithm in our results on sine-Gordon and A(/> 4 one-loop kink 
mass shifts. More importantly, the asymptotic approach has been applied to other systems supporting 
kinks considered in the Literature with a very high level of confidence and effectiveness. These results 
are even more remarkable because, to the best of our knowledge, there are no other ways of computing, 
even approximately, the one- loop kink mass shifts in these models. As a final remark on the advances 
reported here we stress a very appealing property of our improved method: very few data are necessary 
for the computation of one-loop kink mass shifts in this way. Only the potential and the vacua connected 
by the kink are needed. The number of terms kept from the series is optimized in parallel to the choice 
of the integration range in the Mellin transform. 

We plan to keep improving the asymptotic method by applying the procedure to more complex 
situations: 
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1. The first target would be the computation of one- loop kink mass shifts when the kinks interpolate 
among vacua belonging to different points in the vacuum moduli space. In the <?5> 6 -model, this is 
always the case and the idea is to incorporate the improvements achieved here to the calculation of 
the (/> 6 -kink mass shift performed in [28] . There are many more models carrying this type of kink 
susceptible to analysis within the framework of the improved asymptotic approach. 

2. In kink-field theories with more than one scalar real field there are lots of very interesting kinks 
that can have one- two- or ^-components different from zero. The one-loop kink mass shifts of the 
celebrated Montonen-Sarker-Trulinger-Bishop model [43]- [44] are computed in [19]. Now we have 
the possibility of remarkably improving the precision of our results. 

Another field theoretical system of this type is a two-component scalar field theory proposed by 
Bazeia et al. in [ ], who identified some isolated kinks. The solitary- wave solutions were unveiled 
by Shifman and Voloshin in [ 16], whereas we related these kinks to critical trajectories in some inte- 
grable or quasi-integrable dynamical systems in [47] . The one-loop kink mass shifts were discussed 
in [33] and again we expect to get closer to the exact answer by applying the improved asymptotic 
approach. 

3. BPS states in the N = 2-supersymmetric (l + l)-dimensional Wess-Zumino model are kink solutions 
of first-order equations derived from holomorphic superpotentials. Derived in [48], they are crucia 
elements in the classification of N = 2-supersymmetric integrable systems. In [49] we discussed 
these special solitary waves from a real analytical point of view in order to compare them with the 
more traditional kinks arising in the models mentioned above. It seems natural to try the one-loop 
computation of mass shifts for these kinks restricting the theory to the bosonic sector: it is known 
that N = 2 supersymmetry forces the quantum corrections to the BPS masses to be zero. 

4. Finally, we look forward, in a longer perspective, to improving the results on one-loop mass shifts of 
self-dual Nielsen-Olesen vortices achieved in [ ] and [ ] in the Abelian Higgs model. Simili modo, 
important improvements in the results on one-loop mass shifts of the planar self-dual topological 
solitons calculated in [52] in the semi-local Abelian Higgs model should be expected. 
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